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It is known that the localization length scaling of noninteracting electrons near the quantum Hall 
plateau transition can be described in a theory of the bosonic density operators, with no reference 
to the underlying fermions. The resulting "Liouvillian" theory has a global supersymmetry 

as well as a hierarchy of geometric conservation laws related to the noncommutative geometry of 
the lowest Landau level (LLL). Approximations to the Liouvillian theory contain quite different 
physics from standard approximations to the underlying fermionic theory. Mean-field and large- 
N generalizations of the Liouvillian are shown to describe problems of noninteracting bosons that 
enlarge the U(l\l) supersymmetry to U(l\l) x SO(N) or f7(l|l) x SU{N). 

These noninteracting bosonic problems are studied numerically for 2 < N < 8 by Monte Carlo 
simulation and compared to the original N = 1 Liouvillian theory. The N > 1 generalizations pre- 
serve the first two of the hierarchy of geometric conservation laws, leading to logarithmic corrections 
at order 1/N to the diffusive large- iV limit, but do not preserve the remaining conservation laws. 
The emergence of nontrivial scaling at the plateau transition, in the Liouvillian approach, is shown 
to depend sensitively on the unusual geometry of Landau levels. 
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The original work of Anderson on localization in noninteracting electronic systems^ proposes a simple test for the 
. existence of extended states in a disordered system. When an electron is added at the origin at t = 0, its mean 
' squared displacement is finite for all times if all states are localized: ((R 2 (t))) < C for some constant C. If there 
, are extended states in the system, then the electronic motion is diffusive for long times, and ((R 2 (t))) ~ Dt. These 
two alternatives are not exhaustive, however. The motion of lowest-Landau-level electrons in a white-noise potential, 
starting at t = from a localized wave packet, satisfies ((R 2 (t))) ~ Ct a for a w 0.79, and numerics suggest that this 
exponent a is a general property of the universality class of the quantum Hall plateau transition. The physical origin 
of this exponent is the "two-parameter scaling" that describes the Hall transition: both the diagonal and transverse 
conductivities are needed to characterize the system^, unlike in the zero-field case. 

Note that the universal quantity a can be defined without any reference to the electron energy or to the electron 
wavefunction. The electron density correlations are already sufficient to obtain a, which is related to the localization 
length exponent v by a = 1 — This insight is the basis of the Liouvillian approach to the plateau transition 
introduced by Sinova, Meden, and Girvin^. An exact theory which contains a can be written down in terms of the 
O ■ bosonic electron density operators, and while this theory has not yet been solved, it has been studied analytically in 
J> , several different limit o 4 i 5 i 6 . Much of the interest in such noninteracting models stems from the fact that the value of 
v from such models is consistent with that found in some experimental samples 7 . Adding short-ranged interactions 
does not modify this valued, but the effects of long-ranged interactions are currently unclear. 

This paper discusses the physical origin of critical scaling and anomalous dimensions in the Liouvillian picture. 
The bosonic density operators p q in the lowest Landau level satisfy a linear first-order evolution equation; averaging 
this evolution equation over disorder cannot be done analytically. However, previous work has considered mean-field 
and large- TV limits where the disorder-averaged theories can be studied analytically (here N is the number of density 
"flavors" , and N = 1 is the physical case) . Part of this paper considers such theories in more detail and shows how 
the large- N theories represent theories of noninteracting bosons, which can be studied using Monte Carlo simulations 
just as for N = 1. 

An important question about the Liouvillian approach is how the seemingly simple Lagrangian of the density 
operators can show the nontrivial a < 1 scaling behavior observed numerically. We argue that long-lived fluctuations 
resulting from a hierarchy of conservation laws in the Liouvillian formulation cause the scaling. This picture is 
supported by analytical and numerical study of mean-field and large- A" approximations which preserve some but not 
all of these conservation laws. Previous work found diffusive behavior (a = 1) for the mean-field and large- N limits, 
with a diverging set of corrections at order 1/N. These logarithmic corrections occur because of the first nontrivial 
conservation law beyond energy and particle number, that of "chop" , a measure of the variation in the particle density. 
Note that these logarithmic corrections occur for a different physical reason than the logarithmic corrections which 
appear in the theory of weak localization 9 . 

There is a hierarchy of higher conservation laws in the N = 1 problem which are not preserved in the N > 1 
generalizations; these higher conservation laws account for the differences wc find numerically between N = 1 and 
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TV > 1. These conservation laws can be expressed compactly in terms of electron translation around polygonal paths. 
They essentially reflect the nonindependence of the density operators in an LLL state: there are many more density 
operators than independent LLL states, so these density operators must be related to each other. (Conservation laws 
in the Liouvillian formalism are also discussed in&, but those explicitly involve the disorder potential and hence differ 
from the geometric conservation laws discussed here.) 

The standard field-theoretic description of the quantum Hall plateau transition is as a nonlinear sigma model with 
topological termili. This description has generated a great deal of subsequent work and clarified how the plateau 
transition fits into the general description of noninteracting localization problems in terms of nonlinear sigma models 
(NLerMs). However, it has so far been difficult to find a clear understanding of how the topological term yields power- 
law scaling in models with noncompact target spaces, and efforts to obtain scaling via e.g. instanton gas calculations 
are difficult to confirm numerically. Supersymmetric nonlinear sigma models into compact spaces in two dimensions 
can be understood quite thoroughlyAi. These difficulties motivate work on other approaches such as the Liouvillian, 
which has already yielded an improved understanding^ of the connection between the quantum Hall transition and 
classical percolation^. 

Now we explain the connection 3 - between the localization length exponent v and the power-law spreading of a wave 
packet discussed above. The critical exponent v describes the divergence of the localization length near the critical 
energy: 



An electron added in a "typical" localized state (not an eigenstate) will project onto states of many different energies. 
The density of states is believed to be nonsingular at the critical energy, so it is reasonable to assume that on average 
the added electron projects equally onto states of different energies. 

The second assumption needed is that electron motion is diffusive in each state up to the localization length at that 
energy, then stops: 



Here D is a diffusion constant with a nonsingular dependence on energy. Then the mean squared displacement, 
averaged over energy, is 



This prediction a = 1 — ^ can be verified numerically, as discussed in section V. 

The many conservation laws for the density operators mentioned above explain several previous results on mean- 
field and large- TV theories, and suggest a picture which is confirmed below by Monte Carlo numerics. A mean-field 
theory gives only diffusive spreading of a wavepacket at long times 3 . A slight modification to this mean-field theory 
allows it to be interpreted as a large-TV limit, where TV is the number of "flavors" of density operators. Then the 
leading corrections in 1/TV include an infinite series of maximally crossed diagrams, whose sum diverges even though 
each individual diagram is finite. The first extra conservation law, of "chop" , explains this divergence in the same 
way that particle number conservation explains the singularity in weak localization. 

Part of this paper shows how the TV > 1 generalizations can be represented and studied numerically as systems 
of noninteracting bosons in random potentials. The density evolution for finite systems for 2 < TV < 8, found via 
Monte Carlo simulations, does interpolate between the analytic result for TV = oo and the physical TV = 1 quantum 
Hall problem. However, higher conservation laws at TV = 1 do not have equivalents for TV > 1, and as a consequence 
the logarithmic corrections at order 1/TV seem unlikely to indicate a power law for N > 1. The numerical density 
evolution for the largest systems studied at TV = 2, which correspond to Liouvillian matrices of size approximately 
10 5 x 10 5 , is found to behave nearly diffusively for long times, suggesting that the true density evolution for TV > 1 
in the thermodynamic limit is diffusive (possibly with logarithmic corrections), rather than scaling with a power law 
slower than diffusion as for the physical case N = 1. 

An important difference between the quantum Hall transition and the superficially similar classical percolation 
transition is that the former seems numerically to have only one nontrivial exponent in the spreading of the particle 
distribution function, while the latter is "multifractal" and has different exponents for different moments of the 
particle distribution function^. Nevertheless, some aspects of the quantum-mechanical spin quantum Hall transition 
are exactly described by sums over classical percolation hullsii One known property of the Liouvillian action (which 
is derived in the following section) is that truncating higher terms in its nonlocal interaction to make it a local 
theory leads to the classical percolation hull problem^: the full quantum Hall problem is essentially the unique 
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((R 2 (t))) = (Dt) 1 --^ <xt a . 
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generalization of this classical problem to a noncommutative space. Most applications of noncommutative geometry 
to the quantum Hall effect have focused on the plateausi^i&, which are currently less opaque than the transitions; 
the integer transition is a major open problem in condensed matter theory which has a simple statement in terms 
of a field theory in noncommutative geometry. In order to be accessible by both condensed matter and high-energy 
physicists, this paper attempts to be self-contained and reviews some previous work in sections II and III. 

The outline of this paper is as follows. Section II reviews the Liouvillian theory set up to calculate the density 
correlation functions. Section III contains analytical results on this theory and its generalizations to multiple flavors 
of density operators. The density operators for a general LLL state are shown in Section IV to satisfy a hierarchy of 
geometric constraints, which in turn implies that the evolution of density operators under the Liouvillian must preserve 
a number of conservation laws. The large- N generalizations of the Liouvillian preserve the first two conservation laws, 
of particle number and chop, but not the higher conservation laws. Section V uses Monte Carlo simulations of finite- 
size realizations of the Liouvillian and its N > 1 generalizations to test the analytic picture from previous sections, 
and ends with a summary of the main conclusions. 



II. REVIEW OF THE LIOUVILLIAN APPROACH 



The electron density of a single-particle state in the lowest Landau level has a simple time evolution which is closed: 
knowledge of the density at one time determines it for all later times. The Liouvillian equation of motion for the 
density operators in a single disorder realization is derived in the first part of this section. Later sections of this 
paper are chiefly concerned with efforts to average this equation over disorder, either numerically or analytically, and 
extract information about universal quantities such as v. 

The noninteracting Hamiltonian of a single 2D electron in a constant magnetic field Bz and random potential v is, 
after projecting to the lowest Landau level (LLL), 



ff = 5>(-q)Ai- (4) 

q 

Throughout this paper, we will assume the strong-field limit and ignore complications from spin and Landau-level 
mixing. The random potential components are independent Gaussian variables: 

«(q) = «(-q)* 
(Mq)» = o 

9-7T7I 2 

««(q)«(q')» = — <W (5) 

The system is on a square of side L with periodic boundary conditions. Each component of q takes on N$ values, 
where = L 2 /(2w£ 2 ) is the number of flux quanta through the torus and I = yJhc/eB is the magnetic length. 
The equation of motion for the LLL-projected density operator /5 q is 

^Pq = [H, Pq\ = ^ «(-<!') [Pq' > Pq] ■ ( 6 ) 

q' 

This commutator turns out to take an especially simple form because of the special geometry of the lowest Landau 
level. 

There is an identity which connects the LLL-projected densities p q to the "magnetic translation" operators Tqtf 

ft, = e-^/V q . (7) 

The magnetic translation operator r q translate a given LLL state by £ 2 q x z. The commutation relation of the r q 
is fixed by the requirement that translation around a closed path pick up an Aharonov-Bohm phase from the flux 
through the path: 



,f 2 pAq, 



r P , r q ] = 2i sin( )r p+q . (8) 



Using this commutation relation and the Hamiltonian J3J, the equation of motion for the r q is 

Tq — i-Cqq'Tq', (9) 
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where the Liouvillian C is 

^ = ^(q-q')sin(^)e-4l q - q 'r (10 ) 

Now this evolution of the density and translation operators can be connected to quantities of direct interest such 
as the disorder-averaged density correlation function. The solution to the equation of motion @ is 

^(*) = E( e " i£ *) qq 'V(o). (11) 

q' 

Let the averaged density correlation function be defined as 

IT(q,t) = -ii^((Tr{p q (t)p_ q (0)})). (12) 

This is the energy integral of the usual energy-resolved density correlation function—. Its significance will become clear 
in a moment. The solution of the Liouvillian equation of motion then gives 

n(q,t)--^ e - £V/2 (((e- iXt ) qq ». (13) 

In order to understand how IT(q, t) is connected to universal quantities such as is, suppose that the electron starts 
at t = in a localized state with given densities /5 q (0) = /o(q)- The function LI(q, t) explains how the initial condition 
/o(q) evolves over time, on average: 

((p q (t)»=n(q,t)/o(q). (14) 

Later this equation will be used to estimate IT(q, t) numerically. Now the density operators near q = can be used 
to extract moments of the density distribution in real space: for example, the squared displacement (R 2 ) — (x 2 + y 2 ) 
at t = is given by, in the continuum limit, 

<*>~(£ + (15) 

For a localized initial state of one electron, /o(q) — > 1 as q — > 0. At long times t, we expect power-law spreading of 
the wave packet: ((R 2 (t))) ~ t 1 ^^ . The corresponding statement in terms of II(q, t) is 

Hence the universal quantity v is contained in the long-time, small-momentum behavior of II(q, t). For simple 
diffusion, we would have had II(q, t) = e~ Dq *, and ((R 2 )) ~ Dt, as expected. 



III. ANALYTIC RESULTS 



The preceding section reviewed how the critical exponent v of the plateau transition, a problem of nonintcracting 
fermions, is contained in the Liouvillian theory of the bosonic density operators. One would like to understand how 
the special features of the plateau transition are reflected in the Liouvillian and cause the bosonic fields to behave in 
such a complex manner. The first part of this section reviews the mean-field and large- N theories of the Liouvillian. 
The supersymmetry technique (reviewed ini£) is used to obtain a field-theoretic representation of the Liouvillian 3 
with £7(1 1 1) supersymmetry. 

The mean-field or large-TV limit of the Liouvillian shows diffusion of the density. Corrections to this limit are 
therefore localizing, i.e., reduce the exponent a in ((R 2 (t)}} - t a to its physical value a » 0.79. The l/N correction 
is examined analytically for evidence of this localizing effect, and the N > 1 theories are studied numerically in the 
following section. As described in more detail in section V, the orthogonal large- TV theory^ corresponds to N flavors of 
bosonic density operators, not to N flavors of (fermionic) electrons; it is hence inequivalent to e.g. the TV-orbital theory 
of Oppermann and Wegner— . The theory also has some unusual features characteristic of noncommutative-geometric 
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theories, e.g., perturbation theory can be organized in such a way that every individual diagram is convergent, without 
any external regularization. 

Some analytical results are obtained for a large- N generalization where the symmetry of the theory is enlarged to 
£7(1 1 1) x SO(N). There are actually many different possible large- N generalizations, some of which are equivalent 
to the original N = 1 theory. We explain in detail the different physical situations corresponding to the orthogonal 
generalization £/(l|l) x SO(N), a unitary generalization J7(l|l) x SU(N), and the most symmetric generalization 
U(N\N), which is actually equivalent to the original N = 1 theory with f/(l|l) symmetry. All these supersymmetric 
theories will become much more concrete in the following section, when they are realized as different conditions on 
the random potentials of a problem of N flavors of noninteracting bosons. 

The second part of this section shows how the quantum Hall dynamics cause the disorder-averaged correlation 
functions to obey a set of geometric identities. The simplest of these is related to the conservation of "chop", the 
variation in the density for states in the LLL. This conservation law is important for understanding how scaling 
appears at the transition beyond mean-field theory. The soft mode which appears at order 1/N in the large- N theory 
can be shown to exist beyond perturbation theory as a result of this conservation law. 

We begin by reviewing the self-consistent Born approximation to H(q,u), which will later be justified as a large- N 
limit. The disorder-averaged density propagator H(q,uj) is represented in perturbation theory to nth order in the 
potential strength v 2 by diagrams with a single density line and n interaction lines. The propagator of the magnetic 
translation operator has dimensions of time and is defined by 

tl(q,Lu) = he 2 e l2q2 / 2 n(q,Lo). (17) 

The Feynman rules for calculation of this quantity follow from averaging the disorder potentials in the Liouvillian 
(|1U|I . The bare electron line is (u) + iS) , and the Liouvillian interaction vertex induced by averaging over disorder 
has an unusual form: 

4nv 2 22, g 2 I 2 
((£p,p+k£ q +k, q » = tfj2 e ~ e k /2 sin(ypAk)sin(-qAk) 

= i . (18) 



p p + k q + k q 

Now we use the above rules to write down the self-consistent Born approximation (SCBA) 3 , which sums all the 
diagrams with no crossing of interaction lines. In a moment we will show how a very similar approximation emerges 
as the exact large- TV limit of the Liouvillian. Using double lines to denote the SCBA propagator, we have 



(all noncrossing) 



q qi q 

(19) 

q q qi q 

The resulting self-consistency equation for the self-energy E(q, to) is£ 

%u) = ?^ E !_!E^(i^. (20) 

V4 ' 1 h 2 L 2 ^ u> + i6-Z{p,u>) V ' 

Numerical solution of (|2UII in the continuum limit N$ — > oo finds a diffusive form for small hwijv and qi: 

1 2 , (21) 

with Dq w 0.965v£/fr. Note that the SCBA self-energy has a different functional form S ~ ivq 2 £/h from that of the 
first term in perturbation theory in v/lo (Si ~ v 2 q 2 /h 2 aj), though both are consistent with the scaling requirement 

E(q 3 w)=w/(«/w,q). (22) 

The above SCBA calculation predicts diffusion at sufficiently long times, or a = 1. We can understand systematically 
what is missing in the above calculation, and hence what causes subdiffusive scaling a < 1, by first finding a limit 
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in which the SCBA result is correct. Suppose that the counting of diagrams is modified in the following way (in a 
moment the supersymmetric formalism will be introduced to find what theory is described by these diagrams). Let 
each density line carry a flavor i = 1, . . . , N, and constrain the flavor indices at each vertex to be paired as: 



—4— —I— (23) 

i j j i 

Simple counting then shows that a noncrossing diagram of n interaction lines with fixed external indices i = 1 has a 
total index degeneracy N n . Any diagram with at least one crossing has at least one constrained choice of index and 
hence degeneracy smaller by order 1/N (actually smaller by order 1/N 2 for this particular theory) in the N — ► oo 
limit. If we redefine v — > v/^/N in order to keep the effective interaction strength finite, then the \axge-N limit of 
this diagrammatic sum is exactly the SCBA result given above. 

The theory which corresponds to this modified diagram expansion has the same ?7(1|1) supersymmetry as the 
original N = 1 theory, plus an additional SU(N) unitary symmetry. It will turn out that a more useful theory for 
understanding scaling at the plateau transition has orthogonal rather than unitary symmetry. Now we review the 
t/(l|l) theory which describes the original N = 1 problem, then consider its different generalizations. 

The density correlation function II(q, oj) obtained from the Liouvillian 3 can be rewritten as the correlator of a 
complex boson field in a field theory with supersymmetryi&. The theory has one bosonic field <p> and one fermionic 
field "0, which appear symmetrically in the action: 



Pi(q,u) = -i j DcjiDcj) J Dip Dip (p q (p q t 

F(u) = -iu J dq (0 q ^ q + V'qV'qH 

/(I, 2, 3, 4) 



-F( u ) 



"qi Vq2 Vq3 ^q4 



1,2,3,4 

2'0q 1 q2 0q3V'q4 + ^qi ^2 ^q 3 V'q* 



(24) 



The effective interaction from disorder averaging is 

/(I, 2,3,4) 



1 

— e 

7T 



qi ~ q4 ' <S(qi 
1 



q2 - q3 - q4) 



sui I -qi Aq 4 I sin I -q 2 A q 3 



(25) 



The effect of the supersymmetry between the bosonic and fermionic fields can be understood in terms of the Feynman 
diagrams defined above: its effect is to cancel all diagrams for II(q, u>) except those which contain a single bosonic 
line running through the diagram. Such supersymmetry is a general feature of the field theories which describe 
noninteracting quantum-mechanical electrons, but the Liouvillian theory 124|) is superficially quite different from 
conventional a- model descriptions, due to the nonlocal interaction (|25p. 

The class of large- N generalizations introduced in ?? are of the form (note that here the definition of Cj is slightly 
modified to reflect standard usage: now Cj 2 rather than Cj is the coefficient of a term in the action) 



F sen (uj) 



F h 



N 

E 



N 

E 

i=l 



dq( 



dqx dq : 



q^q 

/(I, 2, 3, 4) 
N 



F ir 



+C2 2 (^ qi ^ 2 q3 q4 + 2Vi qi ^ 2 < 3 < 4 + <^ q2 «) 
+C3 2 (^ qi < 2 « + 2«<< + K^LKJ 



\ 



C2 



\ 



c 3 2 ^ V- (26) 

% j j i i i j j i j i_ j 

The physical content of these generalizations will become clearer in Section V, when we discuss the noninteracting 

boson problems that generate them upon disorder averaging. In order to reduce to the original theory at N = 1, 
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we impose ci 2 + C2 2 + C3 2 = 1. The properties of the generalizations are reviewed only briefly here, since more 
details are in£. With C3 2 = f,ci 2 = C2 2 = 0, the only diagrams which survive have just one electron flavor running 
all the way through. In this case the theory for any N is just equivalent to the N — 1 theory, and has a full 
U(N\N) supersymmetry, as each index can be rotated independently (this equivalence for the fully symmetric theory 
is discussed more generally in 11 ). 

Other values of the c, give theories which for finite N are physically different from the N — 1 theory, and are 
analytically solvable in the large- N limit. The mean-field theory of ?? corresponds to the large- N limit of the ci 2 = 1 
theory, which breaks the U(N\N) symmetry down to SU(N) x f/(l|l). This theory was discussed earlier in this 
section. The ci 2 = C2 2 — 1/2 theory has a logarithmically divergent 1/N correction^, whose physical explanation is 
discussed in the following section. For this theory the symmetry is further broken to SO(N) x {7(1 11). The N — > 00 
limit is still diffusive, although the diffusivity is reduced by a factor \J~2 from the unitary case. 

The physical meaning of the N > 1 generalizations is that they correspond to problems of N different flavors of 
bosons moving in random potentials. The bosons are noninteracting and boson number is conserved, as expected since 
the theories remain supersymmetric. The U(N\N) case corresponds to random potentials which never scatter bosons 
of one flavor into bosons of another flavor: then these are trivially equivalent to the C/(l|l) case since the flavors are 
independent. The difference between the orthogonal SO(N) x {7(1|1) and unitary SU(N) x f7(l|l) generalizations, 
for example, is in the correlations between the random potentials. The numerical results in Section V confirm that 
these theories are nontrivial generalizations with physics that interpolates between the N = 1 and N — 00 cases. 

IV. CONSERVATION LAWS AND POLYGON IDENTITIES 

The density operators, which are the fundamental quantities in the Liouvillian approach, must satisfy a very large 
number of conservation laws. If there are states in the LLL, there are Ni different density operators, which 
cannot be assigned values independently: only some values of the density operators correspond to actual LLL states. 
A simple example is that there is no one-electron LLL state with uniform density (i.e., (tplpoip) = 1, (ip\pq\ip) — for 
q ^ 0). The constraints on allowed values of density operators will turn out to be important in understanding how 
approximations to the exact Liouvillian theory. 

A convenient way to represent the constraints on density operators is in terms of geometric "polygon identities." 
Essentially the Liouvillian time evolution must conserve a huge number of quantities, for any realization of disorder, 
in order that the density operators always correspond to an allowed LLL state. The first two conservation laws are 
relatively simple and are preserved by the mean-field and large-TV theories, but higher conservation laws contain 
geometric phase factors and are not always conserved. 

The first conserved quantity is simply the particle number: Co = po = 1- The second is the "chop" C\. for all LLL 
states, 

Ci = X>|T q |$Mr_ q |$ - J2 MtMI 2 = N*- (27) 
q q 

Here the sum is extended over the iV 2 , different momenta in the LLL. This conservation law means that all single- 
particle LLL states have a certain amount of variation of the density, in addition to having fixed total density. It 
prohibits the uniform density state mentioned above, since then C\ = 1. Note that the expectation values make 
this different from the trivial statement ^ r q r_ q = ^2 = These conservation laws clearly have consequences 
for correlation functions of the Liouvillian theory, as explained in more detail below. Now it is shown that Co and 
C\ are just the first two of a hierarchy of independent constraints on the density operators; each constraint induces 
additional restrictions on correlation functions of the Liouvillian. 

Consider translating an electron around a triangle (Fig. ^) by successive action of r qi , r q2 , and r_ qi _ q2 . In doing 
so the electron picks up a phase factor set by the signed area of the triangle, 

exp(i$) = exp(i£ 2 qi A q 2 /2). (28) 

Note that C\ can be thought of as summing over 2-gons, where the electron is first translated by a vector £ 2 q x z and 
then in reverse. Such a 2-gon encloses no area, so $ = 0. Similarly Co can be thought of as summing over the 1-gon 
of zero momentum. The generalization of Co and C\ to triangles is 

C 2 ee £ <r qi )(r qi )(r_ qi _ q2 ) e -^^^/ 2 = TV 2 . (29) 
qi,q2 

Here all three expectation values are taken in one LLL state. Note that periodic or antiperiodic boundary conditions 
must be imposed consistently in this sum when qi + q 2 lies outside the original set of momenta, in order to 
preserve the correct commutation relations. 
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FIG. 1: The conserved quantity C2 is denned as a sum over triangles such as the above, weighted by the phase factor e l *. 

A quick proof that these quantities take fixed values for all states in the LLL starts from the lemma that for any 
state \tf>), the sum over its N4, 2 translations, properly normalized, forms a resolution of the identity: 



E 



1. (30) 



This can be shown explicitly using e.g. the Landau basis of LLL states, in which 

X>|7ql7><e|r- q |/3) = N^ a p5 ie . (31) 
q 

Then the commutation relations and (|30|1 can be used to prove by induction the polygon identity 



C„^ V TT(r qi > \{T_r Je- i *=m. (32) 




Here $ is the signed flux through the area of the polygon formed by the q^. 

These conservation laws imply that the iV-body density correlation functions have singularities at zero total momen- 
tum and small frequency. Conservation of particle number Co implies, of course, that the one-body density correlation 
introduced previously has a singularity at zero momentum: 

ft(q = 0,w) = — J— . (33) 
to + 10 

Conservation of the chop C\ implies a singularity in the two-body correlation function at zero total momentum, for 
any initial momentum, but summed over final momentum: 

V ft 2 (q, -q, qi , - Ql , w) = —J— . (34) 
t-^ to + 10 

qi 

Similar singularities exist in the three-body and higher correlation functions because of the higher conservation laws 
mentioned above, although now the geometric factor e** must be included in the sum over correlation functions. 

The singularity of the two-body correlation function C\ is responsible for the logarithmic corrections found at order 
l/N, where N is the number of density flavors, in the following way. First, the singularity becomes a diffusive pole 
for small total momentum. This is the diffusive pole that appears in the diverging ladder sum which enters into the 
maximally crossed diagrams for the one-body correlation function. These diagrams appear at order l/N, where N 
is the number of density flavors, causing the logarithimic corrections discussed mm and the preceding section. The 
conservation of C\ even in the large- N generalizations (see below) explains why there is a diffusive pole from ladders 
with two density lines, corresponding to jour fermionic lines rather than the two lines in the diagrams leading to the 
weak-localization singularity. Hence the logarithmic divergence at order 1 /N has a different physical origin than the 
logarithmic divergence in weak localization. 

An important difference between Co and Ci on one hand, and the higher conservation laws on the other, is that 
the N > 1 generalizations still preserve analogues of Co and Ci : 
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C l = E< r qX r -q> ( 35 ) 

i,q 

are all conserved, where i = 1,...,N. There is no such generalization for the quantities C„, n > 2. The next 
section shows how numerical simulations can be used to study the N > 1 generalizations in the same way as the well 
understood N = 1 case. The numerical simulations suggest that preserving the entire hierarchy of conservation laws 
C n is important for obtaining the true anomalous (i.e., not diffusive) scaling. 



V. NUMERICAL RESULTS 

Section ITTI derived the Liouvillian equation of motion for the magnetic translation operators r q , 

T q (t) = -^£ qq <T q ^). ( 36 ) 

q' 

This evolution is in the Heisenberg representation: i|36|) also describes the evolution of the expectation values (T q (t)) 
in a specific state, from the initial values (T q (0)). 

Recall that the Liouvillian C is an N<j, x N^ 2 matrix for a torus of N$ states in the LLL, while the Hamiltonian 
is only of size N$ x Na. Hence direct numerical calculation of the Liouvillian equation is grossly inefficient 
compared to evolution of the ordinary Schrodinger equation —ihtp = Hip. Previous numerical studies of the Liouvil- 
lian approach&SiSI have therefore used the Schrodinger equation to evolve ip(t): then the expectation values of the 
translation operators, and hence densities, are just 

(r q (f)) - m)K\m)- (37) 

So the evolution of the N^ 2 x A^-dimensional Liouvillian problem can be found exactly from the reduced Hamiltonian 
problem of size x N<j>. The Liouvillian L is therefore quite redundant: for example, its N^ 2 eigenvalues are just 
all the differences of pairs of eigenvalues of H . 

Now let us consider a generalization of the Liouvillian equation of motion to N flavors of translation operators 
r q , and associated density operators p q = e~ q 1 ^ 4 r q : 

q' 

In a moment we will define the generalized Liouvillian C 1 ^ . Before there were N^ 2 different operators T q which had to 
be evolved; now there are NN^, 2 operators. This simple generalization to N flavors will turn out to have far-reaching 
consequences for the behavior of (|38fl . The key difference is that for TV > 1 there is no underlying Hamiltonian of 
which C tJ is the Liouvillian; this explains how the N > 1 generalizations of the Liouvillian differ from previous work 
on large- N generalized Hamiltonians. 

The generalized Liouvillian C lJ for N > 1 differs from the original C only in having multiple disorder potentials v 13 : 

£^ = |^ ( q-qO S m(^M)e^^'l 2 . (39) 

There are several possible constraints on the potentials v %1 . For instance if v lJ = for i ^ j, then "flavor" is conserved 
and clearly the dynamics are the same as for N = 1. We will concentrate on the case of independent Gaussian random 
variables: 

»«(q) = v ij {-q)* = v ji (q) 

«««(q)» = 

2-7T?; 2 

«^(q)w W (q')» = -^S^faSn + SuSjk). (40) 

This will turn out to yield the orthogonal generalizations mentioned in Section III. The unitary generalizations can 
be obtained by a weaker constraint on the random potentials: instead of v lJ (q) = v v (— q)*, one imposes only 
v ij (q) =v ji (-q)*. 

With this choice of potentials, the equations of motion l|38fl for N > 1 do not have a reduced representation in the 
same way as the N = 1 problem: the generalized Liouvillian for N > 1 is not in general the Liouvillian matrix of 
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any smaller matrix. However, note that (|38|l is still the equation of motion of a noninteracting disordered problem, 
and that the Hilbert space is only polynomial in the system size, rather than exponential as for an interacting 
problem. Another important observation is that the N > 1 problems (|38(l . when disorder-averaged, become exactly 
the supersymmetric SO(N) xU(l\l) theories discussed in Ref£ and section III. So the meaning of the supersymmetry 
of those theories is now clear: they describe A flavors of noninteracting real bosons, which for N = 1 are density 
operators of noninteracting fermions. For A > 1, it is shown below that the bosons cannot be interpreted as density 
operators of some underlying fermions. 

Now we return to the original N = 1 problem and introduce the numerical methods which will later be applied for 
A > 1. The system is on a torus pierced by A^ flux quanta: the Landau basis for the N<j, states of the LLL is 

v rn 

Here -L/2 < x,y < L/2. For definiteness fix A^ odd: then N$ = 2Q + 1 and k = —Q,...,Q in (gTJ. There 
are Nj, independent values of momentum in one direction, spaced by ^ = \J 2i^ j I 2 N ^ iq x = n x y / 2ir/N ( f ) , with 
n x = -Q,...,0,...,Q. 

Now imagine starting the system at t = in the state t/^q, which is localized along the line x = 0. Over time, the 
wavefunction will spread out under the influence of the random potential; this spreading can be studied numerically 
to estimate a = 1 — (2z/) .. The Hamiltonian for a given random potential u(q) is 

H = J2 K-q)P q = £ t,(-q)r q e-^ 2 / 4 . (42) 
q q 

The matrix elements of the translation operators in the basis Ij41ll arc 6 

O'K.nJfc) = e^ n *n y+ 2^ s{j _ k _ % mQd N ^ (43) 

The approach used here starts with a localized state and evolves it under the Hamiltonian (|42() to calculate the mean 
squared displacement over time. This approach is almost as efficient as direct diagonalization of the Hamiltonian for 
the single-flavor problem N — 1, but its real numerical advantage is for the large- A generalizations introduced below, 
where direct diagonalization is unavailable. 

To calculate displacements, we use the expression for {x 2 } in the LLL for a torus with N$ flux quanta: 

{x y- 2 = ^ + £ e- 2 ^(-D"Ko + r_„, )_ 
6 ^— ' 7rn^ 

n=l 

This form may be regarded as a discrete approximation to the second derivative with regard to q x which appears in 
the continuum limit (|15|) . For the state ipo localized around x = 0, we have t Kj o = 1 and (a; 2 ) — > £ 2 /2 as N$ — » oo. 
The rotationally invariant state of lowest angular momentum tp r as —> oo has (x 2 ) = (y 2 ) = I 2 . 
We define a quantity D cst (t) which estimates the diffusion constant over the interval (t/2,t): 

((x 2 (t))) {{x 2 (t/2))) 
D cs t(t) = • (45) 

The motivation for using the interval (t/2, t) is to weaken the effect of the ordered initial condition at t = 0, while 
still averaging over a long interval in time to decrease the statistical fluctuations. 

Results for D est (t) for N = 1 are shown for system sizes from = 81 to = 441 in Fig. El After a short initial 
period, D cst is observed to decline as a power-law until the density begins to reach the boundary of the system: 

D est (t) - i-°- 21±0 01 a w 0.79 ± 0.01. (46) 

This estimate of a corresponds to 2.27 < v < 2.5, consistent with the accepted value 21 . Even for fairly small system 
sizes N > 200 there is a significant region where power-law scaling is observed. Results on the N — 1 system can also 
be used to estimate the scale of finite-size effects for N > 1. 

The main lesson of the above is that it is possible to estimate a from quite small systems even without diagonalizing 
the Hamiltonian completely. It is sufficient to be able to observe the evolution of one initial state for many different 
disorder realizations. Such evolution is numerically fast even for the N > 1 problem, where direct diagonalization 
of the sparse A^ 2 x A^ 2 Liouvillian matrices would require a very large amount of memory. The results in |3 are 
similar to those reported for one system size in&, but differ in that we have used the full Hamiltonian (|42|l without 
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FIG. 2: The decrease of the effective diffusion strength D cst {t) with time for three different system sizes. Time is measured in 
dimensionless units vt/hi, and the diffusion strength is normalized to the mean-field value. The solid line shows D cx t~ ' 21 , 
the value quoted in the text. 

truncation, and used D cst l|45|l to measure displacements. There is also a recent detailed study of the displacement 
scaling in long-range correlated potentials^. 

Now the N > 1 generalizations (j38|l can be considered using the same technique. The initial values of the r 
operators are 

<,(()) = Mm • (47) 
Then after time t, the disorder-averaged correlation function n(q,£) is is estimated, following 114(1 . as 

n(27rfcx/L,t)H(^ }>- (48) 

In order to calculate a, which is determined by the second derivative of II(q, i) at q = 0, we recall that the formula 
for ((x 2 )) on the torus l|44(l is just an approximation to this derivative in a bounded system. That is, a is estimated 
through 

•^' / " < '- i »-'y + ^"»~^ (49) 

6 '-^ Trn z 

n—l 

For large N$, s(0) = i for the initial condition (|47|l . Note that it no longer makes sense to speak of a wavefunction 
ip(t) and displacement (ip\x 2 \ip}. The densities T q (t) cannot generally be obtained as expectation values in any state 
r(t) via (r(t)\T q \^(t)). 

In order to compare N = 2 results with N = 1 results, the disorder strength v of the N — 2 problem must be scaled 
by l/y/2 relative to the strength of the N — l problem. This is an example of the 1/y/N rescaling required to obtain 
a finite theory as N — ► oo in section III. Now the generalized Liouvillian equations of motion (|38|) are integrated 
numerically, and the conservation laws of section IV can be used used to check the accuracy of this integration. 

Another check on the N = 2 numerics is the following test, which also clarifies the difference between the N = 2 
and N — l problems. There are three independent disorder potentials at each momentum q: i> n (q), u 12 (q), and 
v 22 (q). With the added constraint v 11 = v 22 , there are "odd" and "even" densities which evolve independently: 
r ± (q) = (Tq 1 ± Tq 2 )/\/2 evolve like the N = 1 problem in the disorder potentials = (v 11 ± v 12 )/V2. Another way 
to express this change is that the theory's symmetry is enlarged from SO(2) x E/(l|l) to E/(l|l) x J7(l|l) by the extra 
constraint. Numerics confirm that this constrained N = 2 theory has D est (t) equal to that of the N — l theory, and 
quite different from that of the true N = 2 theory for sufficiently long times. 

We now summarize briefly the main conclusions. The large- N generalizations of the Liouvillian theory of the 
quantum Hall effect are supersymmetric theories that describe problems of bosons moving in a disorder potential. 
They preserve some but not all of an infinite hierarchy of conserved quantities in the Liouvillian theory, reflecting the 
interdependence of the momentum operators. The N > 1 generalizations interpolate between the N—l quantum Hall 
case and the mean-field theory at N = oo, and can be studied numerically using Monte Carlo simulations. We find 
that an initially localized wave packet spreads diffusively for sufficiently long times, although the diffusion constant 
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FIG. 3: The evolution of {x 2 ) in units of £ 2 for a finite torus system of = 49 flux quanta, for different values of the flavor 
index N. The N — oo result is obtained from analytic solution of the mean- field equations; the N = 1 solution from direct 
evolution under the Hamiltonian; and N = 2, 3, 4, 8 from direct evolution under the generalized Liouvillian. 
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FIG. 4: The evolution of D ea t(t) in units of Dn=oo for a finite torus system of N$ = 225 flux quanta, for N = 1 and N = 2. 
Note that while the N = 1 diffusion constant continues to decrease, the N — 2 diffusion constant initially decreases but at long 
times begins to level off. 



is significantly reduced from the mean-field value. This study confirms that large- N generalization of the Liouvillian 
formalism is a controlled expansion to calculate some properties of the problem, but suggests that the full hierarchy 
of conserved quantities is required to obtain the complete quantum Hall plateau transition. 

The spreading of an initially localized state with time is shown for various N in Fig. |3J As expected, the finite- N 
theories interpolate between the physical quantum Hall case (N = 1) and the mean-field theory at N = oo. For this 
set of numerics, the system size (49 flux quanta) is too small in order to determine whether the N > 1 theories are 
diffusive (a = 1) or show anomalous scaling: the wave packet reaches the system boundary before the asymptotic 
scaling regime is reached. 

For the N = 2 case, system sizes up to 225 flux quanta can be considered in order to study whether the evolution 
of a wave packet is diffusive at long times. Fig. ?? shows the timc-dcpcndcnt diffusion strength D{t) extracted from 
simulations for N = 1 and 225 states in the LLL, compared to D(t) for N = 2 for the same system size (now the 
Liouvillian matrix is of dimension 2(225) 2 ). Note that for the longest accessible times, the N = 2 problem appears to 
show D(t) leveling off, so (i? 2 ) oc t, rather than D(t) falling to zero as a power-law. We note in passing that the time 
evolution of a state under such large Liouvillian matrices can be obtained relatively rapidly because the Liouvillian 
matrices are quite sparse; finding all eigenvectors and eigenvalues, which is quite simple for the N — 1 Hamiltonian of 
225 states, would require a great deal of computer memory and time for the Liouvillian of the corresponding N = 2 
problem. 

Note added: A preprint by V. Oganesyan, J. Chalker, and S. Sondhi ( cond-mat/0212232) posted soon after this 
manuscript's submission contains a similar conclusion about the behavior of the generalized Liouvillian, reached using 
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different techniques. 
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